Available online at www.sciencedirect.com 
ScienceDirect 


Journal of Power Sources 163 (2007) 679-687 


JOURNAL OF 


www.elsevier.com/locate /jpowsour 


Anode activation polarization on Pt(hk l) electrodes 
in dilute sulphuric acid electrolyte 


R.F. Mann *, J.C. Amphlett, B.A. Peppley, C.P. Thurgood 
Department of Chemistry and Chemical Engineering, Royal Military College of Canada, 
P.O. Box 17000, Station Forces, Kingston, Ontario K7K 7B4, Canada 


Received 30 June 2006; received in revised form 2 October 2006; accepted 3 October 2006 
Available online 17 November 2006 


Abstract 


Proton exchange membrane (PEM) fuel cells have been under development for many years and appear to be the potential solution for many 
electricity supply applications. Modelling and computer simulation of PEM fuel cells have been equally active areas of work as a means of 
developing better understanding of cell and stack operation, facilitating design improvements and supporting system simulation studies. The 
prediction of activation polarization in our previous PEM modelling work, as in most PEM models, concentrated on the cathode losses. Anode 
losses are commonly much smaller and tend to be ignored compared to cathode losses. Further development of the anode activation polarization 
term is being undertaken in order to broaden the application and usefulness of PEM models in general. 

Previously published work on the kinetics of the hydrogen oxidation reaction using Pt(ik l) electrodes in dilute HSO; has been examined and 
further developed for eventual application to the modelling of PEM fuel cells. New correlations for the exchange current density are developed for 
Pt(1 00), Pt(1 10) and Pt(1 1 1) electrodes. Predictive equations for the anode activation polarization are also proposed. In addition, terminology 
has been modified to make the correlation approach and, eventually, the modelling method more easily understood and used by those without an 
extensive background in electrochemistry. 

Crown Copyright © 2006 Published by Elsevier B.V. All rights reserved. 
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1. Introduction 


The background for the present paper has been recently 
presented [1]. Our long-term aim is to further develop our gener- 
alized steady-state electrochemical model, the “GSSEM’, which 
consists of a simple, one-dimensional, model of a PEM fuel cell 
(PEMFC) based on 


V = E + acta + Nact,c + Nohmic + Neonc,a + Necone,c (1) 


Accepting the fact that PEMFC are, in general, not isothermal 
and also typically experience spatial variations in current den- 
sity, two-dimensional and three-dimensional models have been 
developed which have provided insight in the resolution of devel- 
opment and design issues. Such models, however, are much more 
complicated and are computationally intensive. Simple models, 
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such as the GSSEM, lend themselves to simple programming 
and the rapid generation of PEMFC performance predictions. 
Our experience, in the decade since the earliest version of the 
GSSEM was published, is that there is still a need for such one- 
dimensional models. 

The specific goal of the present paper is to further evalu- 
ate a major body of previously published work which studied 
the hydrogen oxidation reaction on Pt(hk/) electrodes in dilute 
H2SO4. This is seen as a useful step in the development of a 
recommended modelling approach for ‘acta, the anode activa- 
tion polarization, in PEMFC. The quantitative prediction of this 
particular loss is generally not given much attention in a PEMFC 
model since it is, normally, much less significant than the cor- 
responding loss at the cathode. There are situations, however, 
where the anode activation polarizations could become signifi- 
cant and a modelling capability is therefore desirable. 

Over the years, there has been an extensive literature on the 
theory and practice of fuel cells in general and proton exchange 
membrane (PEM) fuel cells in particular. In addition, there is 


0378-7753/$ — see front matter. Crown Copyright © 2006 Published by Elsevier B.V. All rights reserved. 


doi: 10.1016/j.jpowsour.2006. 10.012 


680 


R.F. Mann et al. / Journal of Power Sources 163 (2007) 679-687 


F 


Nb 


rds 


RDE 


Nomenclature 

b Tafel slope, 2.3RT/aF (V dec!) 

B Levich constant for RDE (A cm mol”! rpm7!" 2) 

B.-V. Butler—-Volmer 

c concentration (mol cm7?) 

CH, concentration of dissolved hydrogen at the reac- 
tion interface (mol Hz cm7?) 

D diffusivity or diffusion coefficient (cm? s7!) 

E total anode polarization, Mact,a + Nconc,a, in Refs. 


[2-7] (V) or thermodynamic cell emf in Eq. (1) 
(V) 


fi(0,00) functions of 6 and 0ọ in the general Butler-Volmer 


equation 
Faraday’s constant (96,487 clbs equ’t7!) 


GSSEM generalized steady-state electrochemical model 


Henry’s Law ‘constant’ (atm cm? mol7!) 
hydrogen evolution reaction (reverse of Eq. (2)) 
hydrogen oxidation reaction (Eq. (2)) 

current density (A cm~?) 

exchange current density (A cm~?) 

number of electrons being transferred for one act 
of the overall reaction 

number of electrons being transferred ‘after’ the 
rds 

number of electrons being transferred ‘before’ the 
rds 

partial pressure of a gas component (atm) 


PEMFC proton exchange membrane fuel cell 


rate-determining step in the reaction sequence at 
the anode 

gas constant (8.314 J mol`! K7!) 

rotating disk electrode 

temperature (K) 


Greek letters 


a transfer coefficient (in the B.—V. equation) 
(defined by Eq. (5)) 

B symmetry factor 

n polarization (i.e. overvoltage or loss) (V) 

0 fractional surface coverage (generally of 
chemisorbed hydrogen atoms) 

v stoichiometric coefficient (number of times that 
the rds step must take place for the overall reaction 
to occur once) in Eq. (5) and kinematic viscos- 
ity (momentum diffusivity) in Eqs. (13) and (14) 
(cm? s7!) 

w rotational frequency of an RDE (rpm) 

Subscripts 

a anode (or, with n, ‘after’ ) 

act activation or actual 

ad adsorbed 

b before 

c cathode 

conc concentration (relating to mass-transfer losses) 


d diffusion-limited (i.e. mass-transfer-limited) 
(often subscript ‘L’ in other literature) 

exp experimental 

fwd associated with the forward reaction 

H hydrogen atom or Heyrovsky reaction 


H2 hydrogen 


H30* hydrated proton (i.e. H*-H20) 

ohmic relating to ohmic (i.e. iR) losses 

pred predicted 

rev associated with the reverse reaction 

T Tafel reaction 

V Volmer reaction 

0 at zero polarization and zero net current condition 


(i.e. at equilibrium) 


a substantial literature relating to the activation polarization of 
the hydrogen oxidation reaction (HOR) for which the overall 
half-cell reaction is 


Hy) > 2Ht + 2e7 (2) 


The reverse of Eq. (2) is the overall half-cell reaction for hydro- 
gen evolution, the HER. There is a very substantial literature 
relating to the HER, some of which has proved useful in devel- 
oping an HOR model. 

The anode reaction occurs in the presence of a catalyst, typi- 
cally polycrystalline Pt if the anode feed is relatively uncontam- 
inated hydrogen (although mixtures such as Pt/Ru may be used 
if there are carbon-containing species in the anode feed). Since 
the surface of a polycrystalline Pt consists of a mixture of vari- 
ous Pt crystal faces (e.g. Pt(1 00), Pt(1 10) and Pt(1 1 1)), greater 
understanding of the behaviour of the HOR on each type of crys- 
tal could be useful in modelling the HOR on polycrystalline Pt. 

The present paper, in Section 2, reviews some literature 
regarding the mechanism and kinetics of the HOR on Pt(1 00), 
Pt(1 10) and Pt(1 11), primarily published [2-7] by what we 
have arbitrarily called the ‘Markovic Group’. Section 3 sum- 
marizes the development of kinetic equations for a number of 
possible ‘ideal’ HOR processes and then applies these equations 
to the analysis of ‘Markovic Group’ results on anode losses for 
the HOR on the various Pt(h k I) electrodes. 


2. Anode activation polarization and the HOR 
2.1. Introduction 


The basic theories of polarizations (or ‘overvoltages’) for 
the HOR and, especially, the HER have been established for 
decades. The following subsections will briefly review some 
previously published work and develop a correlation approach 
for the anode activation polarization that is a compromise of 
being as mechanistic as possible, yet as simple to understand and 
apply as possible. The goal is to develop a relationship linking 
the anode activation polarization to the current density and the 
other operating parameters of the fuel cell. 
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2.2. The Butler-Volmer equations 


As previously summarized [1], the common starting point 
for the development of an expression relating the activation 
polarization, act or simply n, and the current density for a mul- 
tistep electrode process involving electron transfer in the rate- 
determining step is some form of the Butler—Volmer equation. 
For example [1], the most general form of the B.—V. equation, 
written for the rate-determining step in the overall reaction, could 
be represented by 


Otwa Fn 


RT 6) 


rey Fi 
| = frev exp |- ma 


i 

io = fiwa exp | RT 
where the parameters ffyq and frey, subscripted for the ‘forward’ 
and ‘reverse’ reactions, are functions mainly of 0 and 6p and, 
therefore, depend primarily on i, T and py,. As the current den- 
sity approaches zero, these f parameters will approach a value 
of unity and the general B.—V. equation simplifies to the more 
common 


=i lex { ewa) exp { Qrey FN \ (4) 
E RT RT 


Similarly, the commonly used empirical parameter, the ‘trans- 
fer coefficient’ œ, has been subscripted for the ‘forward’ and 
‘reverse’ reactions. Since both the HER and HOR consist of 
various combinations of single-step reactions, various authors 
(such as Bockris and Reddy [8] and Gileadi et al. [9]) suggest 
that a be defined as follows for such multistep processes: 


= p(n — np — Na) + nb 
v 


(5) 


Definitions of symbols and subscripts are in the Nomenclature. 
The HOR is normally based on the consumption of one molecule 
of hydrogen so that ‘n=? applies. 

The symmetry factor, symbol 6, is commonly used to repre- 
sent the ratio of ‘distance across double layer to summit’ to 
‘distance across whole double layer’ (requiring, therefore, a 
value between 0 and 1). While 6 is commonly assumed to have 
a value of 0.5, a may or may not have that value. 


2.3. Butler-Volmer approximations 


2.3.1. The ‘high-polarization’ B.—V. approximation: the 
Tafel equation 

This is based on the forward reaction being much greater than 
the reverse reaction so that Eq. (3) simplifies to 


H i ( i ) es j ( i ) 

Maca = © | — | In -] = o - 

one aF fiwaio twa Ftwaio 
(6) 


This is the Tafel equation and, if nacta is plotted against 
log(i/ftwa), [2.303RT/afwaF] is the “Tafel slope’, typically sym- 
bol b. 


2.3.2. The ‘low-polarization’ B.—V. approximation 
In certain special cases, further simplifications are possible. 
If fewa and frey are ~1 (i.e. as expected in the micropolarization 


region), if the sum of twa and Oey is 1 (i.e. if Grey is equal to 
{1 — afwa}) and if awa is ~0.5, the general B.—V. equation, Eq. 
(3), simplifies to [1]: 


F 
i= ip H (7) 


2.4. The hydrogen oxidation reaction (HOR) 


2.4.1. Fundamental steps in the HOR 

Various possible steps in the overall HOR have been sug- 
gested. Using ‘M’ to represent an active Pt atom and ‘e~’ to 
represent an electron, these steps are: 


Gi) Ho +MsM-H2 
(ii) M-H2 +M & 2M-H 
‘Tafel 


with (i)+(ii) being the often 


Hy +2M & Haa + Haa. 


step’, simply 


Gii) H20 +M & M-H20 
(iv) M-H2 +M-H20 & M.H +M +e` +H30O* [ion-atom for- 
mation] 


with (i) + (iii) + (iv) being the ‘Heyrovsky step’, often simply 
H? © Ht+Hyagte-. 


(v) M-H+H,0 } H30*+M+e— 
(vi) H30* & H* + H20 [the dehydration step] 


with (v) + (vi), the ‘reverse discharge Volmer step’, often simply 
Hag @ H+ +e. 

It is important to note that the Tafel reaction and the Hey- 
rovsky reaction are in parallel and either one or both of these 
reactions may occur. The Volmer reaction, on the other hand, 
must occur. 


2.4.2. Kinetic expressions for the individual Tafel, 
Heyrovsky and Volmer reactions 

Kinetic expressions (in terms of current density 7) for the 
main reaction possibilities summarized above are as follows, 
essentially using the nomenclature in Breiter’s Eqs. (5)—(7) [10]: 


e For the Tafel reaction: 


ip = kr(l — Oy)’ an, — k-th (8) 
e For the Heyrovsky reaction: 
: annk 
= kyl — 0 
ig H( HAH, exp | RT 
a_-HnFk 
—k_ — 9 
Hay+ Oy exp | RT | (9) 
e For the Volmer reaction: 
. ayn avn 
iy = kyOn exp | n | k-v(1 — Oy aye exp |- = | 
(10) 
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Depending on the catalyst used and the reaction conditions, 
different combinations of parallel and series reactions are pos- 
sible as well as different rate-determining steps. The value of 
the Tafel slope, b, can, in some situations, gives some indication 
of the mechanism followed. Considering the range of literature 
values [3] for the HOR Tafel slope on Pt(h k 1) (e.g. 0.028, 0.037, 
0.074 and 0.112 V dec™! at 274K anda pm, Of 1 atm), a variety 
of mechanisms are apparently operative on different Pt catalysts. 


2.4.3. Kinetic expressions for various combinations of the 
Tafel, Heyrovsky and Volmer reactions 

A general HOR reaction scheme could be called “Tafel (T) 
and/or Heyrovsky (H) (in parallel) followed by Volmer (V)’ (in 
series). Arbitrarily, the T and H reactions could each be at one 
of four relative rates: ‘fast’ (£), ‘slow’ (s), ‘very slow’ (vs) or 
‘zero’ (0). A large number of combinations of {T and/or H} + V 
reactions would, therefore, be theoretically possible. Some of 
these combinations could lead to relatively simple rate equations 
while more complex combinations could lead to very complex 
rate equations. 

The simplest ‘ideal’ HOR processes are “Tafel followed by 
Volmer’ (i.e. zero Heyrovsky) and ‘Heyrovsky followed by 
Volmer’ (i.e. zero Tafel). Either the first or the second reac- 
tion could, ideally, be the rate-determining step (rds) so that the 
derivation of complete rate equations could, therefore, lead to 
four possible ‘ideal’ mechanisms: 


Mechanism I: ‘rds’ Tafel followed by ‘fast’ Volmer; 
Mechanism II: ‘fast’ Tafel followed by ‘rds’ Volmer; 
Mechanism III: ‘rds’ Heyrovsky followed by ‘fast’ Volmer; 
Mechanism IV: ‘fast’ Heyrovsky followed by ‘rds’ Volmer. 


In general, an expression for i(@) is derived from equations 
for the ‘slow’ (rds) reaction and 69 is determined by applying 
an equilibrium analysis to the ‘fast’ reaction. 


3. Evaluation of Butler-Volmer parameters from 
published ‘Markovic Group’ anode polarization data 
for the HOR on Pt(hk D) 


3.1. Introduction to the ‘Markovic Group’ results 


A number of publications [2—7] reported on a series of single- 
crystal RDE studies on Pt(1 00), Pt(1 1 0) and Pt(1 1 1) that were 
carried out by what we have arbitrarily called the “Markovic 
Group’. This set of publications appeared to contain internally 
consistent measurements of Cy, io and anode polarization over 
a wide range in current density, at several temperatures from 1 
to 60°C, in 0.05 M H2S0; electrolyte and for the three Pt(h k 1) 
single-crystal electrodes. The data were reported as ‘adjusted 
to a hydrogen gas pressure of | atm’, the situation that we are 
calling ‘standard conditions’. 

This ‘Markovic Group’ body of work has been evaluated in 
detail, particularly the results at the lowest temperature studied 
(274 K), and the results are summarized in the following subsec- 
tions. In some cases, our proposed parameter values are close to 
those reported by the original authors but in several cases con- 


siderably different parameter values and different conclusions 
are now being proposed. 

The published E values in the Markovic Group publications 
[2-7] represented the total anode polarization, nota, essentially 
the sum of Nact and Neonc. It was decided to extract the E(i) data 
sets from the figures in their publications and to convert these 
to diffusion-corrected Nact,a(i) data sets. Detailed analyses were 
then carried out on these nact,a(i) data sets. 


3.2. Extraction of E(i) data sets 


The published polarization data of the Markovic Group 
appeared in several forms. Using one of their references [3] 
as an example, their Fig. 3 presented ‘raw’ ‘i versus E’ RDE 
plots at 274 K, their Fig. 4 presented ‘E versus log(mass-transfer- 
corrected current densities)’ at 274 K and their Fig. 5 presented ‘i 
versus F’ plots for the micropolarization regions of the 274, 303 
and 333 K data sets. These figures were enlarged and internally 
consistent E(i) data sets were extracted for all three catalysts, 
over the complete range of current density for 274 K and over 
the low current density, micropolarization, range for 274, 303 
and 333 K. The necessary values of limiting current density, their 
ig, were calculated using their Eqs. (1) and (2), the details of this 
calculation being in Section 3.3. 


3.3. Correction of E(i) data sets to Nact(i) data sets 


As described above, this correction involved the application 
of the correlations recommended by the Markovic Group authors 
for quantifying Nconc,a and subtracting it from their E, total anode 
polarization, values to get a new set of Nact,a Values. Their results 
were all obtained with an RDE and they proposed the following 
correlations [3]: 


RT i 
Neonc,a = Nd = OF In | 1 — a (11) 
where 
ig = Bey, o? (12) 


For use in Eq. (12), they presented values of Bcy,, units 
(A cm7? rpm™!?), as follows: 


274K 5.6 x 1075 (experimental) 
298K 6.54 x 1075 (theoretical) 
333K 8.8 x 1075 (experimental) 


The agreement between theoretical and experimental values 
was considered sufficiently good by the authors so that the appli- 
cation of Eq. (11) was straightforward. It was necessary in our 
work, however, to also estimate a value of the Bcy, parameter at 
303 K and, for this, the variation of Bcy, with temperature was 
investigated further. 

A number of things were considered. With emphasis being 
put on the theoretical prediction of B, the Levich constant, a 
starting point was the following [11]: 


B = 0.62nFD*/~!/6 (13) 
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Nact,a (V) 


log i (A/cm?) 


Fig. 1. Anode activation polarization as a function of current density at 274 K 
(data derived from the Markovic Group E(i) data [3]). 


Eq. (13), after the incorporation of the conversion factor 
“1 rpm = 30/z rad’ [2], becomes 


B = 3.88 x 104D? v16 (14) 


Eq. (14), with consideration of their proposed 25°C values of 
D and v [2], temperature coefficients for D and v from literature 
values for water, their proposed experimental values of Bco (i.e. 
our Bcy,) at 274 and 333 K, and evaluation of their cg via our 
proposed Henry’s Law constants for 0.05 M H2S0;4 electrolyte 
[12], becomes 


>| (15) 


B=7.2 x 10° exp -> 
T 


The resulting Nact,a(d) data sets, derived from their 900 rpm data, 
are summarized in Fig. 1 for the complete experimental range 
of current density and Table 1 for the micropolarization range 
of current density. It turned out that the Ncone component of their 
experimental E ranged from about 20 to 50% at this rotation 


rate even in the low current density region. Fig. 1 and Table 1 
data were then used to evaluate ig and the other Butler-Volmer 
parameters as described in the following sections. 


3.4. General comments re-evaluation of Butler-Volmer 
parameters from the ‘adjusted to py, of 1 atm’ results of the 
Markovic Group 


(a) In general, all of the kinetic equations describing Mech- 
anisms I-IV include terms such as 6/09, (1 — 9)/(1 — 90), 
PH>/b PH, and ay+ /pay+, initially unknowns. 

(b) For estimation of ip in the Tafel region, an expression for fry 
as a function of current density must be developed before 
Eq. (6) can be applied. 

(c) For estimation of ig from the micropolarization region 
data, summarized in Table 1 and discussed in Section 2.2 
during the introduction of the ffyq and frev parameters, a 
‘limit’ approach was initially taken (when i — 0, 0 — 6 and 
ay+ — payqt+) so that all of the above ratio terms approach 
a value of unity at i=0. The ‘limit’ forms of the various 
equations were then used to estimate ig(Z) from the (i) data 
sets, the sets of ig(Z) values then being plotted against i and 
extrapolated to zero i to obtain the true ‘ip observed - 

(d) With ig known, expressions for ffwq and frey were then devel- 
oped empirically to fit all the n(i) data. 

(e) It is common in the literature to express the derived values 
of ig as a function of temperature on an Arrhenius plot. 
In general, however, iọ depends on a number of variables 
(e.g. rate ‘constants’, adsorption ‘constants’ and cy, ), all of 
which have their own dependence on temperature, so that 
the overall effect on an Arrhenius plot may not be a simple 
linear ‘In ig versus 1/7” plot. 

(£) The Markovic Group authors also implied that the same a 
value applied at all the temperatures studied, indicating that 
essentially the same mechanism must apply from 274 to 
333 K. 


Table 1 

Pt(hk 1) nacta(i) data sets estimated from the E(i) micropolarization data of the Markovic Group [3] 

Pt(hkl) 274K 303K 333K 
ix 103 (A cm~?) n x 103 (V) ix 103 (A cm~?) n x 103 (V) ix 103 (A cm~?) nx 103 (V) 

Pt(1 10) 0.0655 0.78 0.362 3.24 0.19 1.575 
0.131 1.54 0.50 4.34 0.414 2.90 
0.263 3.00 0.724 6.03 0.647 4.07 
0.491 5.90 1.00 7.95 0.900 5.00 

Pt(1 00) 0.05 1.30 0.241 3.44 0.333 3.07 
0.1 2.60 0.483 7.10 0.632 6.08 
0.2 5.30 0.655 10.24 0.912 9.50 
0.3 8.00 0.776 12.60 1.053 12.7 

Pt(1 11) 0.05 3.10 0.095 1.91 0.12 1.92 
0.1 6.20 0.19 3.79 0.24 3.82 
0.15 9.30 0.362 7.59 0.466 7.61 
0.2 12.6 0.517 11.4 0.69 11.3 
0.25 15.6 0.638 15.4 0.836 15.4 
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(g) As will be seen in the following sections, conventional 
Butler—Volmer equations apply in the Mechanism II and 
Mechanism III analyses of the Pt(1 1 1) and Pt(1 00) results. 
In Mechanism I, the rds does not involve the release of an 
electron so that a ‘conventional’ Butler—-Volmer equation 
does not result in the analysis of the Pt(1 1 0) results. 


3.5. Pt(1 00) results from Markovic Group 


3.5.1. Introduction 

Their analysis [3] of the E(i) results at 274K led to two 
Tafel regions, an intermediate region with a Tafel slope of 
0.037 and a proposed ‘2.3RT/3F’ expression for b plus a 
high current density region (E>0.04 V) with a Tafel slope of 
0.112 and a proposed 2(2.3RT/F) expression for b. This led to 
their conclusion that the process was ‘rds Heyrovsky plus fast 
Volmer’ (i.e. our Mechanism IID) “for the case of small coverage 
by Haq at low overpotentials and ‘maximum’ surface cover- 
age by Haq at high overpotentials”. From their mass-transfer- 
corrected plot of the 274K high current density data, an io 
of about 1.1 x 107° resulted. They then used the simplified 
B.-V. expression, Eq. (7), for the analysis of their microp- 
olarization results, reporting the following exchange current 
densities: 


i9,274K = 0.36 x 10-3 
i0,303K = 0.60 x 1073 
i0,333K = 0.76 x 10-3 


Our analysis of the derived (i) data sets, as summarized in the 
next section, led to somewhat different conclusions. 


3.5.2. Mechanism and analysis of the 274 K Tafel region 
results for Pt(1 00) 

We also concluded that this process was a Mechanism III 
application, rds Heyrovsky plus fast Volmer, but further deter- 
mined that it could be described over the complete range of 
reported current density by the following, derived from Eqs. (9) 
and (10): 


ouo [GE (i) a} 
vet 1-60) \ bpn RT 
0 ayt a_HnF 
- (a) (aorta 
60/7 \ bapt RT 


The ‘2i9 p’ term may also be referred to as ‘ig observed’: 
In the Tafel region (i.e. at higher current densities), the above 
would simplify to 


i X 2io4 = PH, exp aE 
0, 1 = Oo bPH2 RT 


where {(1 — @)/(1 — 90)}(pu, /b Pn, ) can be represented by fiwa- 
The most useful Tafel plot is, therefore, 7 (i.e. Macta) versus 
log(i/ftwa). The parameter ffyq was evaluated by an iterative 
correlation of our Fig. 1, 274K, Pt(1 00) data at high current 
density until the same value of ig was obtained for all the ‘high 
i? data sets. The results, as plotted in Fig. 2, were an ig of 


(16) 


a7) 


0.1 


0.08 


-3.5 -3.3 -3.1 -2.9 -2.7 -2.5 -2.3 -2.1 -1.9 
log[i/fiwa] (A/cm?) 


Fig. 2. Modified Tafel plots at 274 K using the refined x-axis parameter: i/ffwa,274 
from Eq. (18) for Pt(1 00) and from Eq. (28) for Pt(1 1 1). 


1.03 x 107? A cm7? and a Tafel slope of 0.109 V dec~!, both 
relatively consistent with the original Markovic Group values 
(as they should be in this case), plus the following approxima- 
tion: 


fwaz © 1 — 0.37 x 103)” (18) 
3.5.3. Mechanism and analysis of the 274, 303 and 333 K 
micropolarization region results for Pt(1 00) 

The general B.—V. equation, Eq. (16), was used as the starting 
point for the analyses of Table | data. From the high current den- 
sity, 274 K, analysis above, @fwa was kept at a value of 0.5 (i.e. 
a= B). With the assumption that the rds for the reverse reaction 
was the same as for the forward reaction, yey becomes equal to 
1.5 G.e. œ= +1). With pp, /ppH, being ~1, frwa is essentially 
(1 — 0)/(1 — 60) and, in Mechanism III, should decrease as cur- 
rent density increases. Similarly, with ay+/pay+ assumed to be 
~l, frev is essentially 6/69 and would be expected to increase 
with increasing current density. Simple expressions for ffwa and 
frev were modified iteratively until consistent values of ig were 
obtained. The results were: 


e For 274K: 
io,obs = 1.03 x 1077 (i.e. iog = 0.51 x 107°) 


fiwa © 1 — 0.370 x 103)” (19) 


frey © 1+ 1.570 x 10°) (20) 


With the values of iọ and ffyq being the same as obtained 
earlier in the Tafel region analysis, it is concluded that the 
Mechanism III assumption is valid for the complete range of 
current density covered by the Markovic Group work. 

e For 303K: 


io,obs = 1.6 x 1073 (i.e. iog = 0.8 x 1077) 
faa © 1 — 0.2(i x 10°) (21) 


frev ~ 1+ 0.4 x 10°) (22) 
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iprea (A/cm? x 10°) 


(0) 0.2 0.4 0.6 0.8 1 1:2 1.4 1.6 
iexp (A/cm? x 10°) 


Fig. 3. Comparison of iprea, current density predicted for Pt(1 00) from Eqs. 
(25) and (28), and iexp from Fig. 1 and Table 1. 


e For 333 K: 
ip.obs = 3-2 x 107° (ie. iog = 1.6 x 1073) 


fiwa © 1 — 0.2(i x 103)” (23) 


fey © 1+0.2(i x 103) (24) 


The original proposed ip values [3], derived from the E(i) data 
and summarized in Section 3.5.1, do not agree with either of 
the above ip values, ig.obs OF io H (=0.5i0,obs), evaluated from 
our derived n(i) data. 


The quality of the new correlations of the Pt(1 00) results is 
demonstrated in Fig. 3 parity plot, where ipreq values, from Eq. 
(16) plus the results tabulated above, are compared to the iexp 
data from Fig. 1 and Table 1. 


3.6. Pt(1 11) results from Markovic Group 


3.6.1. Introduction 

Their analysis [3] of their Pt(1 1 1) E() results at 274 K led to 
a Tafel slope of 0.074 V dec™!, implying an unusual œ value of 
0.735, a value not representative of any of the ‘simple’ mecha- 
nisms. They concluded that the process involved a parallel ‘fast’ 
Tafel and ‘fast’ Heyrovsky feeding Haq to an rds Volmer, i.e. a 
mixed Mechanisms II + IV. Their mass-transfer-corrected Tafel 
plot indicated an ig of about 0.2 x 107°? Acm7?. They then used 
the simplified B.—V. expression, Eq. (7), for the analysis of their 
micropolarization results, reporting the following exchange cur- 
rent densities: 


i9.274K = 0.21 x 1073 
i9,303K = 0.45 x 10-3 
i0,333K = 0.83 x 1073 


Their i9.274K values from the Tafel analysis and the micropo- 
larization analysis essentially agreed, lending support to their 
correlation approach, but the œ value of 0.735 was a complica- 
tion to further analysis of the process. 

Our analysis of the derived n(i) data sets, as summarized in 
the next section, led to somewhat different conclusions. 


3.6.2. Mechanism and analysis of the 274 K Tafel region 
results for Pt(1 I 1) 

We concluded, unlike the original authors, that this process 
was a Mechanism II application, fast Tafel plus rds Volmer, and 
could be described over the complete range of reported current 
density by the following, derived from Eqs. (8) and (10): 


: : 0 Atwan F 

l = 10,V bo exp RT 
1-0 aH+ = _ Grey F (25) 
1— Oo baH+ P RT 


Now 6/69 can be represented by ffwa and {(1 — 0)/(1 — 4%)} 
(a+ /bap+) can be represented by frev. 

For this situation, 09 is determined by the Tafel reaction equi- 
librium between H2 and Hag. This ‘io,y’ may also be referred to 
as ‘IO observed in a Mechanism II analysis. 

In the Tafel region the above would then simplify to 


a F 
i © iov fv exp {ae \ (26) 


The Tafel approximation for the activation polarization, there- 
fore, is 


oar (220887) m E ia es 


Afw F io 


so that the most useful Tafel plot is nact versus log[i(ffwa)]. The 
fiwa term was estimated by an iterative correlation of our Fig. 1, 
274 K, Pt(1 1 1) data at high current density until the same value 
of iọ was obtained for all the ‘high 7’ data sets. The result, as 
plotted in Fig. 2, was an ig value of 0.25 x 107°? A cm~? and 
a Tafel slope of 0.109 Vdec™!, both differing from the origi- 
nal Markovic Group values. This Tafel slope indicates an awa 
value of 0.5, consistent with a Mechanism II reaction. The above 
results were obtained with the following approximation: 


a 


fewa.274 © 1 — 0.14(i x 10° (28) 


3.6.3. Mechanism and analysis of the 274, 303 and 333 K 


micropolarization region results for Pt(1 1 1) 
The general Mechanism II B.—V. equation, Eq. (25), was used 


as the starting point for the analyses of Table | data. From the 
274 K high current density analysis above, afwq was kept at a 
value of 0.5 (i.e. œ = 6). With the assumption that the rds for the 
reverse reaction was still the Volmer reaction, an Qyey value of 
0.5 (i.e. œ= p) is again obtained from Eq. (5). The parameter 
fiwa is essentially 6/69 and, in Mechanism II, should decrease as 
current density increases. Similarly, with ay+/pa,y+ assumed to 
be ~1, frey is essentially (1 — 0)/(1 — 89) and would be expected 
to increase with increasing current density. Simple expressions 
for ffwa and frey were modified iteratively until consistent values 
of iọ were obtained. The results were: 


e For 274K: 


io,obs = 0.38 x 1077(iLe. ioy = 0.38 x 107°) 
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Fig. 4. Comparison of iprea, current density predicted for Pt(1 1 1) from Eqs. 
(25) and (28), and iexp from Fig. 1 and Table 1. 


fewa as Eq. (28) (as developed from the Tafel region analysis), 
fiwa and frey both ~1 in the micropolarization region, Eq. 
(28) giving values of ffwa just below 1 and further correlation 
giving frey values averaging just above 1. 

With the values of iọ being nearly the same for both the 
Tafel region and the micropolarization region analyses, it is 
concluded that the Mechanism II assumption is valid for the 
complete range of current density covered by the Markovic 
Group work. 

For 303 K: 


io,obs = 1.32 x 107-3 (ie. ioy = 1.32 x 1073) 


fwa © 1 — 0.01 x 103)” (29) 
frev © 1+0.05(i x 10°) (30) 
e For 333 K: 


id,obs = 1.82 x 107°G.e. ioy = 1.82 x 1077) 


a 


fewa © 1 — 0.04(i x 10° (31) 


frey © 1+ 0.0125 x 10%) (32) 


The above expressions for ffwa and frey do all give values 
very close to | at these low current densities but, with the 
forward and reverse reaction rates being quite comparable in 
magnitude, the small decreases in ffyq and the small increases 
in frev as the current density increases do significantly improve 
the validity of the ig evaluation. 


Our values of iọ proposed above are significantly different 
than those proposed in the original publication [3] as are our con- 
clusions regarding the mechanism for the HOR over Pt(1 1 1). 

The quality of the new correlations of the Pt(1 1 1) results is 
demonstrated in Fig. 4 parity plot, where ipreq values, from Eq. 
(25) plus the results tabulated above, are compared to the iexp 
data from Fig. 1 and Table 1. 


3.7. Pt(1 10) results from Markovic Group 


3.7.1. Introduction 

Their analysis [3] of their Pt(1 10) E@) results at 274 K led 
to a Tafel slope of 0.028 V dec™!, implying an « value of 1.94. 
This value was close to the ‘2?’ expected in a Mechanism I, ‘rds 
Tafel and fast Volmer’, process and they concluded that this was 
the correct Pt(1 10) mechanism. Their mass-transfer-corrected 
Tafel plot indicated an ig of about 0.46 x 107°? Acm™?. They 
then used the simplified B.—V. expression, Eq. (7), for the anal- 
ysis of their micropolarization results, reporting the following 
exchange current densities: 


i0,274K = 0.65 x 10-3 
i0,303K = 0.98 x 10-7 
i9,333K = 1.35 x 1073 


3.7.2. Mechanism and analysis of the 274 K Tafel region 
results for Pt(1 10) 

We concluded, as did the original authors [3], that this was 
a Mechanism I situation except that that Eq. (33), derived from 
Eqs. (8) and (10), is our recommended starting point: 


mafia] 


This is of the general form of a Butler—Volmer equation and can 
be rearranged into the form of a Tafel equation, as follows: 


RT | i l 1 — 6p V 
n= In|1—-- 

2F 10,T 1-0 
A Tafel line with slope b = 2.303RT/2F should be expected, indi- 
cating an apparent œ value of 2. At temperatures in the range 
274-296 K, this would therefore give values of b in the range 
0.0272-0.0294 V dec™!, close to the 0.028-0.031 values often 
reported for ‘room temperature’ results for the HOR on Pt. 

Similar to what was done in the Pt(1 0.0) and Pt(1 1 1) sections 

earlier, the parameter {(1 — 6)/(1 — 6o)}? will be represented by 
the symbol fi ; 9 so that the modified Tafel-like equation is 


BI 1 i i | (35) 
n= n 7 
2F fiı1ioio,T 


The 274K n(i) data in Fig. 1 were correlated with an ig of 
1.0 x 107° and an fi 10 given by the following empirical expres- 
sion: 


(34) 


1.05 
fiio.274K = 1 + 0.62(i x 10°) (36) 


3.7.3. Mechanism and analysis of the 274, 303 and 333 K 
micropolarization region results for Pt(1 1 0) 

As i— 0, 0 —> 0o and fi 19 > 1 so that Eq. (33) simplifies to 
z i 

1 — exp{—2F7/RT} 


i0,T (37) 


Eq. (37), with the micropolarization data in Table 1, can be used 
to estimate ig 7, compensating for any variation of ‘apparent’ ip 7 
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Fig. 5. Comparison of iprea, current density predicted for Pt(1 10) from Eqs. 
(33) and (36), and lexp from Fig. | and Table 1. 


with i by extrapolating to io,r at zero i. The following results, in 
Acm~?, were obtained from this analysis: 


io,T,274K = 1.0 x 107° 
10,T,303K = 1.35 x 1073 
i0,T,333K = 1.50 x 107° 


These values of io,r are from 54% to 11% higher than the values 
proposed by the Markovic Group [3] and are slightly non-linear 
on an Arrhenius plot. 

These io,r values were then substituted back into Eq. (35) to 
permit the determination of the following empirical fi 1 o expres- 
sions in addition to the still-valid Eq. (36): 


fi10,303K = 1+ 0.6247 x 10°) (38) 
fi10333K = 1 + 1.16(i x 10°) (39) 


The application of Eqs. (33), (36), (38) and (39), in conjunction 
with the ig values proposed above, in the prediction of current 
density is demonstrated in Fig. 5 parity plot. 


4. Summary and conclusions 


It is recommended that the most general form of kinetic equa- 
tion for each process, represented by Eqs. (16), (25) and (33), 
be used as the starting point for the analysis of polarization 
data for the hydrogen oxidation reaction over Pt(1 00), Pt(1 1 1) 
and Pt(1 10). In general, the functions such as (1 — 0)/(1 — 60) 
and 6/09 are not equal to unity so that they should not be 
dropped from the general kinetic expressions. In general, our 
Eq. (7), which was used by the original authors to evaluate io, 
is only appropriate for data correlation in very special circum- 
stances. 

We concluded, like the original authors, that the high current 
density HOR process on Pt(1 0.0) was a Mechanism III applica- 


tion, ‘rds Heyrovsky plus fast Volmer’, but, unlike the original 
authors, we further determined that the diffusion-corrected n(i) 
data in Fig. 1 and Table 1 could be correlated over the com- 
plete range of reported current density by the Mechanism III 
equations, Eqs. (16)—(24). 

We concluded, unlike the original authors, that the HOR on 
Pt(1 1 1) appears to be well represented, over the complete cur- 
rent density range, by a ‘Mechanism II’ process, “fast Tafel fol- 
lowed by rds Volmer’ and the experimental diffusion-corrected 
n(i) data are well correlated by Eqs. (25) and (28)-(32). 

The HOR on Pt(1 10) appears to be well represented by a 
“Mechanism I’ process, ‘rds Tafel followed by fast Volmer’ and 
the experimental diffusion-corrected (i) data are well correlated 
by Eqs. (33), (35), (36), (38) and (39). 

The values of ig now being proposed average over 50% higher 
than the values originally proposed by the Markovic Group 
since, in all cases, the removal of the neonc,a component from 
their anode polarization, Æ, led to a lower value of nact,a and, 
therefore, a higher value of io. 

In the correlation of the Pt(1 00) and Pt(1 11) results, the 
various œ values were equal to either 6 or to 6+1 and, in all 
cases, a B value of 0.5 appeared appropriate, consistent with the 
common recommendation in the literature. 
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